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Abstract. 

The gonihedric Ising Hamiltonians defined in three and higher dimensions by 
Savvidy and Wegner provide an extensive, and little explored, catalogue of spin models 
on (hyper)cubic lattices with many interesting features. In three dimensions the k = 
gonihedric Ising model on a cubic lattice has been shown to possess a degenerate low- 
temperature phase and a first order phase transition, as well as interesting dynamical 
properties. The dual Hamiltonian to this may be written as an anisotropic Ashkin- 
Teller model and also has a degenerate low-temperature phase as a result of similar 
symmetries to the original plaquette action. 

It is possible to write an alternative dual formulation which utilizes three flavours of 
spins, rather than the two of the Ashkin- Teller model. This still possesses anisotropic 
couplings, but all the interaction terms are now four spin couplings and it acquires 
an additional, local gauge symmetry. We investigate this alternative dual Hamiltonian 
using zero temperature and mean-field methods together with Monte-Carlo simulations 
and discuss its properties and the relation to the Ashkin- Teller variant. 
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1. Introduction, and One Dual 

The gonihedric Ising models are a family of lattice spin models defined by Savvidy 
and Wegner whose spin cluster boundaries are weighted to mimic the worldsheets of 
a gonihedric string action which was also originally formulated by Savvidy jlj. If the 
string worldsheet is discretized using triangulations, this action may be written as 

S = ^\jti-X j \0(a ij ), (1) 

(ij) 

where 6(a.ij) = \ir — QZij\, a%j is the dihedral angle between the neighbouring triangles 
with a common edge (ij) and |JQ — Xj\ are the lengths of the embedded triangle edges. 
The word gonihedric is a neologism resulting from the combination of the Greek words 
gonia for angle (referring to the dihedral angle) and hedra for base or face (referring to 
the adjacent triangles). 

The idea of using generalized Ising models to investigate a gas of random surfaces 
was employed by Karowski |2] and Huse and Leibler [3] amongst others and was 
described in some detail by Cappi et.al. in [I]. The plaquettes of the geometrical 
spin cluster boundaries give the random surfaces of interest and one can relate an Ising 
Hamiltonian with nearest neighbour (ij), next to nearest neighbour ((ij)) and plaquette 
[ijkl] terms 

PH = - Ji ^ ~ J z ° ia S ~ Js a i a i a k°h (2) 

{ij) ((ij)) [ijkl] 

whose partition function is given by 

z = j2^M-m (3) 

M 

to a surface partition function for a gas of surfaces 

Z = ex V (-(3 A A(S) - f3 I I(S) - (3 C C(S)) (4) 
{s} 

where S is a configuration of spin clusters, A(S) is the total area of boundary plaquettes 
in spin clusters, I(S) is the number of edges shared by four plaquettes and C(S) is the 
number of contiguous plaquettes meeting and right angles (i.e. "edges"). The two sets 
of couplings are related by 

(3 A = 2 Ji + 8 J 2 
/3 C = 2 J 3 - 2 J 2 

fr = - 4J 2 - 4J 3 (5) 

which can be seen by enumerating the possible spin values around the local plaquette 
configurations that give rise to edges, intersections and boundaries. 

The gonihedric surface and Ising models are characterized by tuning the area 
coupling Pa to zero, which fundamentally changes the critical characteristics of the 
spin Hamiltonians. This fixes the ratio of nearest neighbour and next to nearest 
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neighbour couplings, so if we set the edge coupling to be (3 and ask that the energy of an 
intersection, given by convention to be (3j+2(3c, is 4/3 k; we arrive at the Hamiltonian (5p] 



H = -2kJ2 a i a 3 + ^Y Gi(T i ~ Y i°i°kOi ( 6 ) 

(v) (m mi] 

where we have extracted an overall factor of for convenience. When k = the 
gonihedric Hamiltonian becomes a purely plaquette term 

H k= q — — o~iO~j(JkO~i (7) 

[ijki] 

where the spins live on vertices and it is the dual(s) of this latter Hamiltonian which 
are the focus of this paper. 

One formulation of the dual to this k = gonihedric Ising model was constructed 
by (Savvidy) 3 in [7], 

H = ~ 2 Y ai °i ~ 2 Y TiTk ~ \ Y W (8) 

(ij) (ik) (jk) 

where the sums are carried out over the orthogonal axes ij, ik and jk on the 3d cubic 
lattice (so each sum is one-dimensional) and the spin variables on the vertices live in the 
fourth order Abelian group. For simulations it is perhaps easier to deal with Ising (±1) 
spins, which may be done choosing rji = <JiJi with (Ji — Ti — ±1, which still satisfies the 
requisite algebraic relations. The net result is that the dual Hamiltonian is given by the 
anisotropic Ashkin- Teller model [8] 

H Ashkin-Teller = ~ ^ Yl GiG j ~ 0~ Ys TiTk ~ 2 Yj (X i <TkT i Tk ' ^ 

(ij) (ik) (jk) 

In the isotropically coupled case the ratio in equ. ^ corresponds to the Z 4 
symmetry point of the standard Ashkin-Teller model which can be seen by rewriting the 
Hamiltonian in terms of the four double spins Si = (±1, ±1) to give the 4-state Potts 
Hamiltonian 

H = -\Yi A6 ^-i) (io) 

(ij) 

where the sum now runs over all three edges orientations. This isotropic Hamiltonian 
has a first order transition in three dimensions |9~1. 



We found in 10 that HAshkin-Teiier, as for its isotropic Ashkin-Teller counterpart, 
displayed a first order transition. Unlike the isotropic model it had a highly degenerate 
ground state (and probably low temperature phase) because of the possibility of flipping 
planes of spins at zero energy cost. Simulations also suggested that its dynamical 
features were rather similar to the original plaquette Hamiltonian, with metastability 
around the phase transition and, possibly, glassy characteristics at lower temperatures. 
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2. More Duals 

Savvidy and Wegner also considered the formulation of higher dimensional equivalents of 



the plaquette gonihedric Hamiltonian and their duals in 11 . Borrowing their notation 
temporarily, a general Hamiltonian in the style of the gonihedric Hamiltonian for d — n 
dimensional hypersurfaces on a d dimensional hypercubic lattice can be written as 

H ai ...a n+1 ( r ) — ~ K ^ ] ^ ] Va 1 ...a k _ 1 a k+1 ...a„ +1 ,a k ( r ) (H) 

ai<...<a„+i,r k 

where the individual terms are given by the products 

Va 1 ...a k ^ 1 a k+1 ...a n+1 ,a k ( r ) — £^«i...a fc _ 1 a fc+1 ...« n+ i ( r ) 

x U ai ...a k _ 1 a k+1 ...a n +i ( r — e a k ) iX^) 

with the LPs being —1 if the dual hyperplaquette belongs to the hypersurface and 
1 otherwise. This Hamiltonian is intended to count the number of hyperplaquettes 
in the d — n dimensional hypersurface which do not continue straight through a 
hyperedge, which is the natural generalization of the three dimensional gonihedric 
plaquette Hamiltonian. 

With this choice of Hamiltonian the energy associated with non-straight hyperedges 
is AK times the length of the bend. To obtain closed surfaces we can again follow the 
three dimensional case and introduce Ising spins a attached to (d — n + 1) dimensional 
hypercubes £l ai ...a n -i- The V terms may then be written as the product over the spins 
at the boundary of two parallel neighbouring (d — n)-dimensional hyperplaquettes 

n 

Vai...an,f}( r ) — J_ J_ CT ai...a k _ 1 a k+1 ...a n ( r ) cr ai...a k _ 1 a k+1 ...a„ ( r ~ e a k ) 
k=l 

x cr a 1 ...a k _ 1 a k+1 ...a„( r ~ &!3)&a 1 ...a k _ 1 a k+1 ...a n ( r ~ &a k ~ 6/3). (13) 

For instance, in four dimensions the resulting Hamiltonian for two dimensional 
hypersurfaces was 

HgLhedric = - ^(^ff) P "Hff)p (14) 

p 

where the sum runs over pairs of parallel plaquettes P in 3d cubes on the Ad lattice and 
the Ising spins a are now located on the centres of the edges in the four dimensional 
lattice. 

Several variants of the duals for such Hamiltonians were discussed by Savvidy and 
Wegner. In the case of two dimensional surfaces in general d the dual Hamiltonian was 
found to be 

H d dual , 2d = ~^^A^r(r>-+e>V(f+e 7 ). (15) 

where there are d(d — l)/2 A spins at each vertex and T spins on each edge. On the 
other hand, the dual Hamiltonian for hypersurfaces of codimension one was given by 

H duai,codimi= ~K* ]\^ a ,^(^^ a Mr + e 1 )A^ ai (r)A^ ai (r + e^) 

a</3,r 7 

(16) 
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with the standard duality relation relating the coupling to the original K 

exp(-2iT) = tanh(2AT) . (17) 

For two dimensional surfaces in three dimensions either formulation may be employed. 
The codimension one variant gives a dual Hamiltonian with three different Ising spins 
{Ai, A 2 , A 3 } at every vertex r, 

Hdiai,codimi = - A a (f)A^(r)A a (f+e 7 )A^(f+e 7 ) (18) 

which at first sight looks rather different to the Ashkin- Teller-like Hamiltonian discussed 
in the introduction which contains only two spins. H^ al codiml can also be derived from 
the two-dimensional surface Hamiltonian of equ. (15) by summing over the T edge spins. 

Reverting to the notation used in the introduction to facilitate comparisons, let us 
take Ai(r) = A 2 (r) = n and A 3 (r) = /ij which allows us to rewrite equ. (18) as 

Hdual = — ^ y %2i <T i <J jl J >il J >j ~ n T i r k^iHh ~ q Vj&kTjTk (19) 

(ij) (ik) (jk) 

where the sums are again carried out over the orthogonal directions ij, ik and jk on 
the 3d cubic lattice and for compactness we have denoted this Hamiltonian as simply 
Hdual- Although the Hamiltonian of equ. (19) contains three spins there is a local gauge 
symmetry 

o-j, Tj, Hi ->■ -a h -n, -Hi (20) 
or alternatively 

a i7 n, ^ ->■ %n, ^ifjii (21) 

where 7j is also an Ising spin, which we would expect to reduce the number of local 
degrees of freedom back to two. We shall see explicitly how this occurs below. 

In the rest of the paper we investigate the behaviour of the dual Hamiltonian Hdual 
using zero temperature and mean field calculations along with Monte-Carlo simulations 
and compare it with both the Ashkin- Teller style dual formulation HAshkin-Teiier of 
equ. ([9]) and the original plaquette Hamiltonian H K=0 of equ. ([7]). 



3. Ground State(s) 

To investigate the ground state (i.e. zero-temperature) structure of the theory whilst 
allowing for possible non-uniform states the Hamiltonian H dU ai may be written as a sum 
over the individual cube Hamiltonians he [4], 

h c = - g ^2 aiCfk ViVi ~ g TiT i ^ k ~~ g °i° k T i Tk ' ^ 

(ij) (ik) (jk) 

where the additional symmetry factor of 1/4 with respect to the full Hamiltonian is a 
consequence of one edge being shared by four cubes. Minimizing the energy of the full 
Hamiltonian can then be achieved by minimizing the cube energy and tiling the full 
lattice with compatible minimum energy configurations of individual cubes. 
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A feature of both H K=0 and its Ashkin- Teller dual HAshkin-Teiier is that it is 
possible to flip planes of spins (which may be intersecting) at zero energy cost, giving 
a peculiar symmetry that lies somewhere between the gauge and the global. For the 
dual Hamiltonian under consideration here the anisotropic couplings mean that it is 
possible to flip planes of pairs of spins at zero energy cost as shown in Fig. ([I]). It is 
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Figure 1. Some possible ground state spin configurations on a cube, the a, t, fi 
values are shown at each vertex. The directions of the anisotropic couplings in the 
Hamiltonian are indicated, as are the faces on which pairs of spins are flipped. 



also possible to flip two or three orthogonal faces on the cube, so tiling the entire lattice 
with such combinations we can see that in addition to the purely ferromagnetic ground 
state we may have arbitrary (and possibly intersecting) flipped planes of pairs of spins. 
In addition, the local gauge symmetry may also be employed to flip all three spins at 
any vertex. 

The ground state structure, and the mechanism of anisotropic couplings which 
allows the plane spin flips, is clearly very similar to that in HAshkin-Teiier 



10 , whose 



possible ground states we show in Fig. ([2]). We can make the relation even more explicit 
by using the gauge symmetry to reverse the signs of the spins on the shaded planes in 
Fig. ([l]b,d) which fixes the third spin fi to be positive in all configurations. The a and 
r spins then have the same configurations as in suitably rotated configurations of the 
Ashkin- Teller dual. 

In summary, the ground state structure of H duai shows an interesting interplay 
between the "semi-local" planar flip symmetry which appears to be a characteristic of 
both the gonihedric models and their duals and a local gauge symmetry. This allows one 
to reduce the effective number of degrees of freedom and recover the ground states of 
HAshkin-Teiier, which not only has anisotropic couplings but a different form of coupling 
(energy-energy rather than spin-spin) in one of the directions. All the couplings in H^uai, 
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Figure 2. Possible ground state spin configurations on a cube for the Ashkin- Teller 
formulation of the dual Hamiltonian. The a, r values are shown at each vertex. The 
directions of the anisotropic couplings in the Hamiltonian are again indicated, as are 
the faces on which spins are flipped. 



on the other hand, are energy-energy but the price to be paid for this is the introduction 
of the new gauge symmetry. 

4. Mean Field 

Applying a mean field approach to systems with non-uniform low temperature phases 
can be done in a similar fashion to the discussion of the ground states by using a 
cube decomposition [21. The total mean field free energy is again written as a sum of 
the elementary cube terms 4>(lc,mc,nc), now given by a sum of energy and entropy 
contributions 

P(f>(lc,m c ,n c ) = - - ^ rnirrijlilk-- ^ n i n kkh 

(ij)cC (ik)cC 
\- 

(jk)cc 

+ 4 + *o ln d + *o + ( x - ^ ln d - 

ice 



+ — 5^(1 + rrii) ln(l + + (1 - m,) ln(l - mj] 

ice 

+ ^ ^[(1 + ln(l + ri,) + (1 - m) In (1 - n<) ] (23) 



ice 
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where the m^n^li are the average magnetizations corresponding to the spins cr^Tj and 
/ij respectively and the In terms are entropic factors. We then minimize the cube free 
energy numerically using the 24 equations 

d(j)(lc,mc,n c ) 



9h (i=i.... 







d(f)(lc,m c ,n c ) = 

d</>{lc,m c ,n c ) = Q , 24 , 

dfli (j=1...8) 

or, more explicitly 

mi = tanh[/3(m 4 / 1 / 4 + m 2 n\ n 2 )} 

m 8 = tanh[/3(m 5 / 5 / 8 + m 7 n 7 n 8 )] 

rii = tanh[/3(n5/i/5 + n 2 mi m 2 )] 

n$ = tanh[/3(n 4 / 4 / 8 + n 7 m 7 m 8 )] (25) 

h = tanh[/3(/ 4 mim 4 + / 5 ri\ n 5 )} 

Is = tanh[/3(Z 4 n 4 n 8 + Z 5 m 5 ms)] 

where, as is our wont, we have labelled the magnetizations on a face of the cube counter- 
clockwise 1...4 and similarly for the opposing face 5 . . . 8, as shown in Fig. (|3|. In 



1 2 




mn 



Figure 3. The labelling of vertices used in writing the mean field equations for the 
cube. The directions of the mean field spin couplings in the Hamiltonian are again 
indicated. 

disordered phases the mean field magnetizations will iterate to zero, whereas in an 
ordered phase they will take non-zero values. 
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As in our previous work with such equations stability issues may be addressed by 
modifying the iterative scheme slightly from 

= ai\m\n k ], mf +1) = 9l [l k , m k , n k ] , nf +1) = h[l k : m\ n k \ , 

(26) 



to 



(fc+i) 
m) 



nf +1) = V ' ,L ' -V: 1 ' "~" (27) 
1 + a 

for a suitable a, to deal with the possibility that it might fail to converge if an eigenvalue 
of dl\ k+r * /dlj , dmf' +V> /dm k or dnf + ^ /dn k is less than — 1 ji . 

The use of coupled equations on a cube accommodates non-uniform solutions in an 
identical manner to the ground state discussion. Any such solution can then be used to 
tile the full lattice. If we solve the coupled mean field equations iteratively a single first- 
order transition is found in the region of /3 — 0.98. The ground state chosen depends on 
the initial seed values for the l,m,n magnetizations. Choosing these to be near +1 will 
lead to the simple ferromagnetic ground state, other values will pick out one of the other 
possible ground states (or ones related to them by the local gauge transformation that 
reverses all of the spin signs at a single vertex). Although the transition temperature 
calculated from the mean field equations for H^uai is different to that observed for the 
mean field equations in the Ashkin- Teller dual ((3 = 0.83) we shall see below that Monte 
Carlo simulations show there is a much closer similarity between the two systems. 





m k , n k ] 


+ alf) 




1 + a 




(9i[l\ 


m k , n k ] 


+ am k J 




1 + a 




(K[l k 


m k , n k ] 


+ an k ) 



5. Monte Carlo 

We carry out Monte-Carlo simulations using 10 3 , 12 3 , 14 3 , 16 3 , 18 3 and 20 3 lattices with 
periodic boundary conditions for the a, r and fi spins at various temperatures with a 
simple Metropolis update. After a suitable number of thermalization sweeps determined 
by the energy autocorrelation time, 10 7 measurement sweeps were carried out at each 
lattice size for each temperature simulated. 

The first observation to make is that there is no signal of a phase transition in the 



magnetizations (a), (r) or (//}, as was also the case with the Ashkin- Teller dual 10 
Unlike the Ashkin- Teller dual, however, the corresponding susceptibilities for H^uai also 
stay constant, whereas in the former they showed a sharp drop at the transition point. 
Given the ability to flip planes of spins and also reverse the sign of all of the spins at a 
vertex in a ground state the absence of a simple magnetic order parameter is no surprise, 
suggesting that these symmetries persist throughout the low temperature phase. 

The phase transition is apparent in measurements of the energy, where there is 
a sharp drop. A plot of the energy is shown for various lattice sizes simulated in 
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P 

Figure 4. The energy for lattices ranging from 10 3 to 20 3 from left to right. The lines 
joining the data points are drawn to guide the eye. Hot starts have been used in the 
various simulations. 



Fig. Q. The energies measured at given inverse temperatures f3 and at given lattice 
sizes are identical to within the error bars for those measured with the Ashkin- Teller 
Hamiltonian. We have already seen that the ground states of H dua i may be mapped to 
those of H Ashkin-Teiier by a suitable choice of gauge transform, but the two Hamiltonians 
may be related more generally by carrying out the gauge transformations 

CTj, T i: ^ ->■ fliCTi, fMiTi, /4 (28) 

i.e. 

0i, T i, Hi -> Vi&i, Hi T i, 1 • (29) 
This relates the partition functions for Hduai and HAshkin-Teiier 
Z= ^ exp[-f3H dua i(a,r,fi)] 

= 2 L * eX P [-pH du al{<T, T,H= 1)] (30) 

W,r} 

= ^ eX P \-^Ashkin-TelleM, t)] . 

W,r} 

Fixing the gauge in Hduai thus corresponds to setting /ij = 1, Vz. The result of carrying 
out this operation, a unitary gauge fixing, is to reduce the partition function for the 
three spin Hamiltonian of Hduai to the partition function for the two spin Ashkin- Teller 
dual Hamiltonian. 

Since Hduai is an "un-gauge-fixed" version of HAshkin-Teiier we would expect the 
sharp drop in the energy seen in Fig. Q to still be indicative of a first order transition 
as it was for HAshkin-Teiier- To confirm this we histogrammed the energy during the 
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simulations, which should display a two-peak structure near the finite size pseudo- 
critical temperature for first order transitions. An example for a 10 3 lattice close to 
its finite size pseudo-critical temperature at (3 C = 1.275 is shown in Fig. We have 
histogrammed the 10 7 measurements of the energy which were carried out after each 
full lattice sweep of the a, r and \i spins. The expected double peak structure is clearly 




Figure 5. The energy histogram from a simulation with 10 7 sweeps on a 10 3 lattice 
near the finite size transition point (3 — 1.275. P(E) is shown on a logarithmic scale. 



visible and is again strikingly similar to that for the Ashkin- Teller dual at the same 
temperature. 

At a first order transition we would also expect to find a non-trivial limit for Binder's 
energy cumulant which is defined as 



1 - 



(31) 



3(£ 2 ) 2 ' 

This approaches 2/3 at a second order transition point and a non-trivial limit at a first 
order point, which is the case here. We can use the finite size scaling of the position 
of this non-trivial minimum to obtain an estimate for C since the expected scaling of 
the pseudo critical value (3 m i n (L) is (3 m i n (L) = (3 C — 0(1/L 3 ) at a first order transition. 
If we plot the measured minima positions for the various lattice sizes against 1/L 3 
we get a reasonable fit to this behaviour with a value of (3 C ~ 1.391(3) and a Xdof 
of 2.25 when the smallest lattice size of 10 3 is dropped from the fits. The estimated 
value of C from this procedure is consistent with the behaviour of the energy jump 
and an estimate from slow cooling in the next section. It is also consistent with the 



estimate obtained using the same procedures for the Ashkin- Teller dual model in 10 of 



1.388(4). These estimates of /3 C are both rather larger than an estimate resulting from 
taking the dual of measurements of the original plaquette action's critical temperature, 
(3 C = 0.54757(63), in 12 . Allowing for various factors of 2 in the coupling definitions, 



an estimate for the dual transition temperature is then given by the standard formula 
as — ln[tanh(/3 c /2)] = 1.32. It is possible that the fixed boundary conditions employed 
in (12), as well as changing the finite size scaling corrections from the standard periodic 
boundary conditions, may have introduced other systematic errors since they have the 
effect of projecting out the ferromagnetic low energy state. The consistency of the 
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1/L 3 

Figure 6. The scaling of the position of the minimum of the Binder energy cumulant 
against the inverse volume. 

estimate for the dual /3 C from the Ashkin- Teller dual with the estimate here for Hduai 
gives some confidence in the quoted values for the dual(s). Taking (3 C = 1.391 gives a 
dual value around 0.51, lower than the direct estimate in [12] but close to an earlier 
estimate in (61. 




6. Dynamics 



One of the most interesting aspects of the original plaquette Hamiltonian is its dynamical 
behaviour. It possesses a region of strong metastability around the first order phase 



transition and displays glassy characteristics at lower temperatures 13, 14 . We 



found that the Ashkin- Teller dual Hamiltonian also appeared to share some of these 
characteristics since it failed to relax to the minimum energy of E = —1.5 when cooled 



quickly from a hot start 10 



Hduai displays identical behaviour under cooling to the Ashkin- Teller dual. We 
consider 20 3 , 60 3 and 80 3 lattices which are first equilibrated at T = 3.0 and then cooled 
at different rates to zero temperature. The energy time series is recorded during this 
process. In Fig. Q we can see that with a slow cooling rate of 5T = 0.00001 per sweep, 
the model still relaxes to a ground state with E = —1.5 for all the lattice sizes. However, 
as can be seen in Fig. ^ with a faster cooling rate of ST = 0.001 per sweep the model 
no longer relaxes to the ground state energy of E = —1.5, but is trapped at a higher 
value, which is around —1.415 for the larger two (60 3 and 80 3 ) lattices. The temperature 
at which the jump in the energy under slow cooling, T ~ 0.72 occurs is consistent with 
the estimate of the transition temperature (3 ~ 1.391 arrived at by other means in the 
previous section. 

The plaquette model also presents very similar dynamical behaviour, which has 



been further investigated in 13 14 , showing a region of metastability around the 
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Figure 7. The time series of energy measurements obtained from cooling 20 3 , 60 3 and 
80 3 lattices from a hot start at T — 3.0 at a rate of ST — 0.00001 per sweep. The 
traces are effectively indistinguishable. 




Figure 8. The time series of energy measurements obtained from cooling 20 3 , 60 3 and 
80 3 lattices from a hot start at a rate of ST = 0.001 per sweep. 



transition point and non-trivial ageing properties at low temperature. Whether this is 



genuinely glassy behaviour or not remains a matter of debate. It was suggested in 14 



that in the plaquette model the supercooled high-temperature ( "liquid" ) phase becomes 
physically irrelevant at a temperature below the observed transition point, giving what 
is effectively an equivalent of the mean field spinodal point. 



Another Dual Gonihedric 3D Ising Model 



14 



7. Discussion 

We have investigated an alternative dual formulation of the plaquette 3d gonihedric 
Ising Hamiltonian of equ. ^ and clarified its relation with the Ashkin- Teller dual of 
equ. ([8]). An additional gauge symmetry in the Hamiltonian Hduai considered here allows 
one to reduce the three local degrees of freedom to two and maps the Hamiltonian onto 
that of the anisotropic Ashkin- Teller dual HAshkin-Teiier- There appears to be no simple 
magnetic order parameter for the low temperature phase in either case because of the 
flip symmetry and, for Hduai-, the additional gauge symmetry. The shared properties of 
the two dual formulations extend to their dynamics, where fast cooling leads in both 
cases to a state which has a higher energy than the ground state(s). 

It is a curious feature of both HAshkin-Teiier and H dua i here that, although they 
are dual to an isotropic model, the anisotropic nature of the couplings plays such 
a fundamental role in determining their properties. We could take "simplifying" 
the Hamiltonian at the expense of introducing further spins and symmetries a stage 
further by disentangling the four spin interactions in Hduai by using an additional 
edge spin a = 1,2,3 acting as an auxiliary field which may be integrated out 
(or, more correctly summed out) to give the codimension one Hamiltonian. This, in 
effect, amounts to using the two-dimensional surface variant of the dual Hamiltonian 



in equ. (15) rather than the codimension one variant that gives Hduai- This gives a 
Hamiltonian of the form 

(ij) (ik) 

-^(^ + ^) ( 32 ) 

W 

which we would expect to present identical behaviour to H Ashkin-Teiier and Hduai- An 
isotropic version of such a Hamiltonian with one flavour of spin would simply be the 
gauge-matter coupling term in the Hamiltonian of the Z2 gauge-Higgs model which 



displays no transition on its own, so the anisotropy of equ. (32) must again play an 
important role in determining the phase structure. 
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